A dynamical model for competing opinions 
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We propose an opinion model based on agents located at the vertices of a regular lattice. Each 
agent has an independent opinion (among an arbitrary, but fixed, number of choices) and its own 
degree of conviction. The latter changes every time it interacts with another agent who has a 
different opinion. The dynamics leads to size distributions of clusters (made up of agents which 
have the same opinion and are located at contiguous spatial positions) which follow a power law, as 
long as the range of the interaction between the agents is not too short, i.e. the system self-organizes 
into a critical state. Short range interactions lead to an exponential cut off in the size distribution 
and to spatial correlations which cause agents which have the same opinion to be closely grouped. 
When the diversity of opinions is restricted to two, non-consensus dynamic is observed, with unequal 
population fractions, whereas consensus is reached if the agents are also allowed to interact with 
those which are located far from them. 

PACS numbers: 05.65.+b,89.75.-k,87.23.Ge 



I. INTRODUCTION 

Statistical mechanics has turned out to be quite suc- 
cessful in modeling many systems whose interaction is, 
in principle, much more complex than those tradition- 
ally studied in physics as, in many cases, the systems 
are made up of agents which are endowed with intelli- 
gence and, therefore, the interaction between them de- 
pends on their decisions Nevertheless, simple sta- 
tistical models have been developed for describing social 
systems drill, economy (BHU? etc. Despite the great com- 
plexity of such real systems, their main properties can be 
reproduced by simple models which retain their underly- 
ing features. 

A great deal of effort has been devoted to developing 
models for describing the properties of systems made up 
of agents with competing opinions p], 0, 0413 • This is of 
great relevance as human conflicts very often arise from 
the simultaneous existence of incompatible opinions in 
populations. Different systems, such as hierarchical soci- 
eties [9] or democratic ones pjj, where the agents follow 
the opinion of the local majority within a group, have, for 
instance, been investigated. Most of these models allow 
the agents to assume only one of two possible opinions. 
Such spin flip models are representative of many real sit- 
uations which offer only two possibilities and, therefore, 
are also of great interest [l2|, ll3| , besides the similari- 
ties with other physical systems. 

The evolution of scientific paradigms has been recently 
modeled in Ref. Q • The slow decline of old ideas and the 
quick adoption of new ones are the main characteristics 
of the model. Those authors find that the dynamics nat- 
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urally leads to the replacement of old concepts by a new 
paradigm, which dominates for a certain period of time, 
until it is gradually replaced by new ideas. 

Inspired by that work, we have developed an opinion 
model where agents are placed at the vertices of a reg- 
ular lattice and interact only with those located within 
a certain range. Each agent has an opinion and its de- 
gree of conviction. In contrast with many other models, 
such as those developed in Refs. 0, HH G3 , for instance, 
the interaction between two agents is strictly local, in 
the sense that it relies only on the agents' properties, 
i.e. their opinions and convictions. Their neighborhood 
has no influence on the interaction. The latter primar- 
ily affects their degree of conviction. More precisely, the 
interacting agents' convictions are affected during the in- 
teraction. If the conviction of one of the agents reaches 
a certain lower bound then its opinion changes to that of 
the opponent. We therefore take into account the diffi- 
culty in persuading someone who has a strong conviction. 
In such a case, it is necessary to change his (her) beliefs 
prior to the acceptance of the new idea. Furthermore, we 
allow the agents to interact with those which are located 
beyond their first neighbors. 

The dynamics of the opinion distribution in popula- 
tions is then studied in the framework of this model and 
the remainder of this work is organized as follows. In 
Sec. [Ill we give a detailed description of the model. The 
results are presented and discussed in Sec. Hill whereas 
concluding remarks are drawn at Sec. IIVI 



II. THE MODEL 



The system is built on a regular mesh of N x hor- 
izontal and N y vertical lines, with periodic boundary 
conditions. One agent is placed at each vertex and 
an opinion C^, among the N Q possible ones, is ran- 



domly assigned to the i-th agent, as well as a positive 
integer which corresponds to its degree of convic- 
tion. An agent located at vertex (fc, I) interacts with 
any of the neighbors located in vertices (k' ,V) ^ (fc,Z), 
where k! = k — r, k — r + l,---,fc + r — 1, fc + r and 
Z' = I — r, Z — r + 1, • • • , I + r — 1, Z + r. The range r is one 
of the model parameters. Thus, each agent has 4r(r + 1) 
neighbors. At each step of the dynamics: 

a) An agent i is sampled with probability proportional 

to d and one of its neighbors is randomly selected 
for interaction, as described below. By doing so 
we assume that the agent's activity is related to its 
conviction. 

b) With probability a, another agent is randomly se- 

lected among the others and its opinion changes to 
any of the N Q possible ones, including its own. This 
procedure represents the replacement of the agent 
by death or substitution due to departure from its 
neighborhood. Its new conviction is selected be- 
tween 1 and the maximum existing value, in order 
not to introduce any bias into the system. 

In step (a) above, nothing is done if the agents have 
the same opinion and one then proceeds to step (b). 

Two agents T and 'j' may interact only if their opin- 
ions are different. In this case, they do it with probability 



p = exp[-A(C min /C max ) 2 ] , 



(1) 



where C m i n (C max ) is the minimum (maximum) between 
Ci and Cj , and A is a parameter which is chosen so as to 
minimize the interaction between agents for which d ~ 
Cj . This is because it is very unlikely that a leader would 
be influenced by another competing leader (two agents 
with large and similar values of C). On the other hand, 
if a leader (large C value) meets an ordinary agent (small 
C value) who has a different opinion, the latter is very 
likely to be convinced by the former. The functional form 
chosen in this work aims at introducing these features 
into the model. If the agents interact, with probability 
d/(d + Cj), the conviction d is increased by one unit 
and Cj decreases by \d — Cj\. Otherwise, Cj increases 
by one unit and Ci decreases by \d — Cj\. If Cj < 
(Ci < 0) then its opinion changes to that of agent i (j) 
and Cj (Ci) is set to unit. 
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FIG. 1: (Color online) Populations of the different opinions 
as a function of time. The different groups have been shifted 
in order to prevent overlap between them. The curves show 



the results for r = 1, A = 1.0, N Q = 5, a = 1.0 x 10" 
N x = N y = 100. For details, see the text. 



and 



obtained using A = 1.0, r = 1, N Q = 5, a = 1.0 x 10 -6 , 
and N x = N y = 100 are shown in Fig. [TJ The pop- 
ulations have been shifted in order to prevent overlap 
between them. For clarity, we have also restricted the 
time scale to 10 6 , in spite of having carried out simula- 
tions up to much larger times, as just mentioned. One 
sees that no opinion dominates the dynamics. They co- 
exist in different proportions and one notices that, very 
often, there is one which is much more popular than the 
others. Its dominance lasts for a relatively short time 
and the popular opinion is replaced by another one. As 
a matter of fact, only a small number of opinions are 
effectively disseminated through the system, the others 
being a small perturbation most of the time. The most 
popular opinions are replaced by the unimportant ones, 
but few opinions dominate the population at the same 
time. 

This conclusion is independent of the number of pos- 
sible opinions N Q as one sees in Fig. [2j which shows the 
number of opinions, for which 



^AWotal > e 



(2) 



where 



III. RESULTS AND DISCUSSION 

At the initialization stage, an opinion 1 < Oi < N Q 
and the conviction degree 1 < Ci < 10 [25] are randomly 
selected and assigned to the i-th agent. The system then 
evolves during, at least, 10 8 steps. A full step corresponds 
to the time interval during which N x x N y intermediate 
steps, as explained in (a) and (b) in SecHIl take place. 

We start by examining the time evolution of the pop- 
ulations of groups with the same opinion. The results 



N 

iVtotal = ( 3 ) 

as a function of time, for e = 0.1 and 0.5. Before carry- 
ing out the sum, the indices i are rearranged so that 
the smallest populations are selected to calculate the 
fractions. The calculations have been carried out for 
N = 50. It is clear from these results that, although 
the composition may change, only very few opinions are 
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actually adopted by the populations. Therefore, the use 
of a fixed number of opinions should not be seen as a lim- 
itation of the model as the system naturally eliminates 
most of the competing opinions and very few of them ef- 
fectively take part into the dynamics. We have checked 
that this conclusion still holds if one uses different pa- 
rameters, such as A = 2.0 or r = 10. 
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FIG. 2: (Color online) Sum of populations, such that 
Ni/Ntot&i > e, e = 0.1 and 0.5. The parameter set is 
the same used in Fig. [TJ except for the number of opinions 
where N Q = 50 is now employed. For details, see the text. 

We now turn to the size distribution of clusters made 
up of neighbors holding the same opinion. The system 
configuration is analyzed at every 10 3 full steps and the 
time average is thus performed. The results are exhib- 
ited in Fig. [3] for different values of the range r. Two 
agents belong to the same cluster if they have the same 
opinion and their grid coordinates obey (fc, I) — (k f , V) = 
(±1, 0), (0, ±1), or (zbl,zbl). The results clearly show 
that the size distribution is very sensitive to the inter- 
action range r. For cluster sizes up to 10% of the total 
system, the distribution becomes steeper as r increases, 
whereas the development of a big cluster, of approxi- 
mately the size of the total system, becomes more and 
more pronounced. Since agents interact only if they do 
not have the same opinion, the borders between clusters 
are the regions of strong activity. For short range interac- 
tions, only agents which are located very close to the bor- 
ders are allowed to interact. Therefore, for small values of 
r, one should expect to observe compact groups of agents, 
who share the same opinion. This should favor the ap- 
pearance of medium size clusters. Indeed, large range 
values would lead to very diffuse borders and, therefore, 
to the disappearance of the coherence among the agents 
which are close to each other. In the limit of very large 
r values, the connected agents would pervade the system 
and the different groups would interpenetrate each other, 
as they would not be compact. This would favor the ap- 
pearance of large clusters during the dynamics, whose 
contribution to the size distribution may also be noticed 
in Fig. 02 Owing to the strict conservation laws, the mul- 



tiplicity of small clusters should then decrease and their 
size distribution would be steeper, as is also seen in Fig. 

El 

This qualitative reasoning is confirmed by the results 
displayed in Fig. [4j where the spatial configuration of the 
clusters is shown at randomly selected moments. Distinct 
opinions are represented by different gray values (color 
online). One sees that, compact groups are indeed formed 
for r = 1, whereas the clusters become more and more 
spatially diffuse as r increases. Thus, our model predicts 
that long range interactions tend to destroy spatial cor- 
relations among opinions, when consensus has not been 
reached and different opinions coexist in the system. 
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FIG. 3: (Color online) Size distribution of clusters made up of 
neighbor agents which have the same opinion. The parameter 
set corresponds to that used in Fig. [2] The power laws are 
best fit to the results, whose exponents j3 — 1.65, 1.90, 2.30, 
and 2.35 are respectively associated with r — 1, 2,5, and 10. 
In the inset, the model predictions for A = 1.0 are compared 
to those for A = 2.0 and r = 1 in both cases. For details, see 
the text. 

The role played by the parameter A is illustrated in 
the inset of Fig. [3j where the cluster size distributions 
obtained with A = 1.0 and A = 2.0, for r = 1, are com- 
pared. The effect on the size distribution is small and is 
more easily noticed at large sizes where one observes a 
slight suppression of big clusters. We have also checked 
that the other observables studied in this work are weakly 
affected if one changes A in the range 1.0 < A < 2.0. For 
the sake of simplicity, we adopt A = 1.0 from here on. 

One may also notice in Fig. [3] that the size distribution 
of clusters, whose size s is smaller than 10% of the system 
size, is fairly accurately approximated by a power law, i.e. 
P(s) oc s~P . The exponent varies with r and corresponds 
to f3 = 1.65, 1.90, 2.30, and 2.35 for r = 1, 2, 5, and 10, 
respectively. We have checked that the asymptotic value 
is reached for 5 < r < 10. 

These results show that the system self organizes, i.e. 
the configurations are reached without the need of the 
external tuning of any parameter, and the power law sug- 
gests the existence of critical behavior. The self organized 



criticality (SOC) has been discussed in different places 
16-20] and has been observed in many different systems 
2CH24j . In our case, this suggests that, on rare occa- 
sions, consensus would spontaneously be reached. This 
state should survive for a while, until conflicting opinions 
are nucleated by the noise described in (b) in Sec. HJand 
the competition among them would restart and another 
opinion would dominate, and so on. 
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FIG. 4: (Color online) Spacial cluster distribution at ran- 
domly selected moments of the dynamics. Agents sharing the 
same opinion are depicted by the same gray values (color). 
The simulation has been carried out using N Q = 5, A = 1.0, 
and N x = N v = 100. For details, see the text. 



In order to investigate whether scale invar iance is, in 
fact, present in the dynamics, we show, in Fig. [5j the 
cluster size distributions for different r values and system 
sizes. The results reveal that, for r = 1, there seems to 
be a characteristic scale, since the largest cluster formed 
during the dynamics does not scale with the system size. 
In fact, as shown by the full circles in the up left frame of 
this figure, the distribution for r = 1 and N x = N y = 200 
is accurately described by P(s) oc exp(— 0.00027s)/s 1,65 , 
which has a clear exponential cut off. These statements 
are weakened for r = 2 and are no longer valid for larger 
range values. More precisely, the power law regime ex- 
tends to larger sizes as the total system size increases 
for r > 2. Actually, there is a range value, between 2 
and 5, for which the power law regime is an adequate 
description of the size distribution, except for very large 
clusters, since finite size effects have to be considered 
for those clusters. Thus, the dynamics leads to two dis- 



tinct scenarios. For small r, compact clusters are formed 
and, occasionally, amalgamate and form a very large one. 
Due to the locality of the coalescence process (for small r 
the interactions take place at the cluster's borders), this 
happens very rarely since it requires strong spatial corre- 
lations. Then, this process leads to the existence of a size 
cut off. On the other hand, when large r values are used, 
the coalescence extends through much larger areas, due 
to the spread of the clusters. It therefore makes it easier 
for correlations to propagate through the entire system. 
Thus, the present model predicts the existence of SOC in 
systems with competing opinions, if the interaction be- 
tween the agents are not restricted to their contiguous 
neighborhood. 
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FIG. 5: (Color online) Same as Fig. [3] for different system 
sizes. The same parameter set has also been employed in 
the present simulations. The full circles in the up left frame 
correspond to P(s) oc exp(— 0.00027s)/s 1-65 . For details, see 
the text. 

Up to this point, we have not investigated the role 
played by the parameter a, which regulates the frequency 
with which an agent randomly changes his opinion. It 
contributes with noise, which prevents the system from 
freezing when consensus is reached. In this sense, the 
model strongly relies on this parameter to ensure an in- 
teresting dynamics. We have found that it also plays a 
very important role in determining the cluster size distri- 
bution. We observed that a = 10 -5 still leads to power 
law regimes for not too large clusters. However, the in- 
variance of the exponent with the system size shown in 
Fig. [5j for r > 2, does not hold in this case. This means 
that, although noise is needed by the dynamics, too much 
noise destroys the scale invar iance, i.e. the agents must 
keep their opinions for, at least, a short while, in order 
to preserve spatial correlations. Since, on the average, 
N x x N y a agents randomly change their opinions at each 
step, one sees that there is no unique value of a that 
would ensure scale invariance for arbitrary system sizes 
since one may always find a size for which too much noise 
is added to the system at each time step, destroying the 
spatial correlations. This shortcoming may be avoided 
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by redefining a as the total rate per step, i.e. propor- 
tional to (NxNy) -1 , so that the desired amount of noise 
is introduced into the system during the dynamics, for 
any system size. 
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FIG. 6: (Color online) Time averaged conviction distribution. 
The parameter set is the same employed in the calculations 
shown in Fig. [3] For details, see the text. 

Since the agents' convictions play an important role 
in the dynamics, as it directly influences their resistance 
to the adoption of new paradigms, we also examine this 
quantity. Thus, Fig. [6] displays the time averaged con- 
viction distribution for different values of r. As in real 
life, most of the agents have a low degree of conviction 
and the system has very few leaders (large C values), i.e. 
the distribution decays exponentially. As expected, for a 
given value of C, the distribution falls off as r increases, 
since longer range interaction allows the agents to en- 
counter others with different opinions more often (as is 
illustrated in Fig. 2J agents with the same opinion tend 
to be closer for small r values and they do not interact). 
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FIG. 7: (Color online) Time evolution of the difference be- 
tween populations with opinion 1 and 2. Except for N — 2, 
the parameter set is the same employed in the calculations 
shown in Fig. [4] For details, see the text. 

We finally examine the dynamics when only two opin- 



ions are allowed. The modeling of such systems is of g reat 
interest and have been extensively studied p], [Io|, H34l3 | 
since there are many situations in real life where binary 
choices have to be made [1]. Furthermore, the above 
results, associated with Fig. [2j also suggest that this sce- 
nario should retain most of the properties of real sys- 
tems. The results of the model simulation obtained for 
N x = N y = 100 are displayed in Fig. [7J which exhibits 
the difference between the populations with opinion 1 
(ni) and opinion 2 (712). As in Refs. [HlGjl, where a non- 
consensus opinion model has been proposed and studied, 
our model allows for the dynamic coexistence of the two 
conflicting opinions with unequal population fractions. 
This is seen in the left panels of this figure which show 
the results for short range interactions between neigh- 
bors, i.e. r = 1 and 2. This is in agreement with Ref. 
[l2j which considers the interaction between the closest 
neighbors. For larger r values, one sees that consensus 
is reached and it lasts for a long period. We have fol- 
lowed the dynamics during much longer time scales and 
confirmed this feature. Owing to the noise introduced 
by the random change of the agents' opinions, regulated 
by the parameter a, the status quo does not last forever 
and, after being the overwhelmingly dominating opinion 
for a long period, its replacement occurs very quickly 
and the other opinion becomes the consensus, and so on. 
Therefore, our model predicts that a transition from non- 
consensus to consensus occurs as the interaction among 
the agents changes from short to long range. 



IV. CONCLUDING REMARKS 

We have developed a model for the dynamics of com- 
peting opinions, which is based on the agents' degree of 
conviction and on the range of the interaction between 
them. It predicts that, even when many different opin- 
ions are allowed, only very few of them are really in use 
by the agents during the dynamics. This is, in fact, ob- 
served in real life when, for instance, at the beginning 
of an election process, many candidates running for a 
political office start with not too different opportunities 
but, after a while, very few dominate the voters' pref- 
erences. The model also predicts that the size distri- 
bution of clusters, made up of agents which are located 
in contiguous spatial positions and share the same opin- 
ion, follows a power law. That distribution is reached 
independently of the initial conditions, i. e. the dynam- 
ics leads to SOC [l6l-[20j, as long as the interaction be- 
tween the agents is not restricted to too close neighbors. 
When only two opinions are allowed, the model leads to 
non-consensus dynamics, which qualitativel y ag rees with 
the non-consensus model proposed in Ref. [12]. On the 
other hand, if the agents also interact with those who 
are located relatively far from them, consensus is quickly 
reached and it lasts for a long time. The dominating 
opinion is occasionally replaced, but there is consensus 
almost all the time. Our model then provides a means to 
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simulate many of the properties of real systems by chang- 
ing a parameter which has a very simple interpretation 
on physical systems, i. e. the range of the interaction 
between the agents. It also contrast with other models 
as the interaction between agents affects their conviction 
in first place and their opinions change only after their 
paradigms have been corroded. 
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